Efficient molecular dynamics using geodesic integration and solvent-solute splitting

We present an approach to Langevin dynamics in the presence of holonomic constraints based on decomposition of the system into components representing geodesic flow, constrained impulse and constrained diffusion. We show that a particular ordering of the components results in an integrator that is an order of magnitude more accurate for configurational averages than existing alternatives. Moreover, by combining the geodesic integration method with a solvent–solute force splitting, we demonstrate that stepsizes of at least 8 fs can be used for solvated biomolecules with high sampling accuracy and without substantially altering diffusion rates, approximately increasing by a factor of two the efficiency of molecular dynamics sampling for such systems. The methods described in this article are easily implemented using the standard apparatus of modern simulation codes

Ergodic properties of quasi-Markovian generalized Langevin equations with configuration dependent noise and non-conservative force

We discuss the ergodic properties of quasi-Markovian stochastic differential equations, providing general conditions that ensure existence and uniqueness of a smooth invariant distribution and exponential convergence of the evolution operator in suitably weighted $L^{\infty}$ spaces, which implies the validity of central limit theorem for the respective solution processes. The main new result is an ergodicity condition for the generalized Langevin equation with configuration-dependent noise and (non-)conservative force

Optic Flow Stimuli in and Near the Visual Field Centre: A Group fMRI Study of Motion Sensitive Regions

Motion stimuli in one visual hemifield activate human primary visual areas of the contralateral side, but suppress activity of the corresponding ipsilateral regions. While hemifield motion is rare in everyday life, motion in both hemifields occurs regularly whenever we move. Consequently, during motion primary visual regions should simultaneously receive excitatory and inhibitory inputs. A comparison of primary and higher visual cortex activations induced by bilateral and unilateral motion stimuli is missing up to now. Many motion studies focused on the MT+ complex in the parieto-occipito-temporal cortex. In single human subjects MT+ has been subdivided in area MT, which was activated by motion stimuli in the contralateral visual field, and area MST, which responded to motion in both the contra- and ipsilateral field. In this study we investigated the cortical activation when excitatory and inhibitory inputs interfere with each other in primary visual regions and we present for the first time group results of the MT+ subregions, allowing for comparisons with the group results of other motion processing studies. Using functional magnetic resonance imaging (fMRI), we investigated whole brain activations in a large group of healthy humans by applying optic flow stimuli in and near the visual field centre and performed a second level analysis. Primary visual areas were activated exclusively by motion in the contralateral field but to our surprise not by central flow fields. Inhibitory inputs to primary visual regions appear to cancel simultaneously occurring excitatory inputs during central flow field stimulation. Within MT+ we identified two subregions. Putative area MST (pMST) was activated by ipsi- and contralateral stimulation and located in the anterior part of MT+. The second subregion was located in the more posterior part of MT+ (putative area MT, pMT)

Order through Disorder: Hyper-Mobile C-Terminal Residues Stabilize the Folded State of a Helical Peptide. A Molecular Dynamics Study

Conventional wisdom has it that the presence of disordered regions in the three-dimensional structures of polypeptides not only does not contribute significantly to the thermodynamic stability of their folded state, but, on the contrary, that the presence of disorder leads to a decrease of the corresponding proteins' stability. We have performed extensive 3.4 µs long folding simulations (in explicit solvent and with full electrostatics) of an undecamer peptide of experimentally known helical structure, both with and without its disordered (four residue long) C-terminal tail. Our simulations clearly indicate that the presence of the apparently disordered (in structural terms) C-terminal tail, increases the thermodynamic stability of the peptide's folded (helical) state. These results show that at least for the case of relatively short peptides, the interplay between thermodynamic stability and the apparent structural stability can be rather subtle, with even disordered regions contributing significantly to the stability of the folded state. Our results have clear implications for the understanding of peptide energetics and the design of foldable peptides

Ecology of the Scorpion, Microtityus jaumei in Sierra de Canasta, Cuba

An assessment of the population dynamics of Microtityus jaumei Armas (Scorpiones: Buthidae) on the slopes south of Sierra de Canasta, Guantánamo Province, Cuba show an increase in activity over the year (≤ 0.05). The activity peak is related to the reproductive period from June to November. The abundance of scorpions was significantly related to density of the canopy and thickness of the substrate

Prefrontal Norepinephrine Determines Attribution of “High” Motivational Salience

Intense motivational salience attribution is considered to have a major role in the development of different psychopathologies. Numerous brain areas are involved in “normal” motivational salience attribution processes; however, it is not clear whether common or different neural mechanisms also underlie intense motivational salience attribution. To elucidate this a brain area and a neural system had to be envisaged that were involved only in motivational salience attribution to highly salient stimuli. Using intracerebral microdialysis, we found that natural stimuli induced an increase in norepinephrine release in the medial prefrontal cortex of mice proportional to their salience, and that selective prefrontal norepinephrine depletion abolished the increase of norepinephrine release in the medial prefrontal cortex induced by exposure to appetitive (palatable food) or aversive (light) stimuli independently of salience. However, selective norepinephrine depletion in the medial prefrontal cortex impaired the place conditioning induced exclusively by highly salient stimuli, thus indicating that prefrontal noradrenergic transmission determines approach or avoidance responses to both reward- and aversion-related natural stimuli only when the salience of the unconditioned natural stimulus is high enough to induce sustained norepinephrine outflow. This affirms that prefrontal noradrenergic transmission determines motivational salience attribution selectively when intense motivational salience is processed, as in conditions that characterize psychopathological outcomes

Ab initio molecular dynamics of liquid water using embedded-fragment second-order many-body perturbation theory towards its accurate property prediction

A direct, simultaneous calculation of properties of a liquid using an ab initio electron-correlated theory has long been unthinkable. Here we present structural, dynamical, and response properties of liquid water calculated by ab initio molecular dynamics using the embedded-fragment spin-component-scaled second-order many-body perturbation method with the aug-cc-pVDZ basis set. This level of theory is chosen as it accurately and inexpensively reproduces the water dimer potential energy surface from the coupled-cluster singles, doubles, and noniterative triples with the augcc-pVQZ basis set, which is nearly exact. The calculated radial distribution function, self-diffusion coefficient, coordinate number, and dipole moment, as well as the infrared and Raman spectra are in excellent agreement with experimental results. The shapes and widths of the OH stretching bands in the infrared and Raman spectra and their isotropic-anisotropic Raman noncoincidence, which reflect the diverse local hydrogen-bond environment, are also reproduced computationally. The simulation also reveals intriguing dynamic features of the environment, which are difficult to probe experimentally, such as a surprisingly large fluctuation in the coordination number and the detailed mechanism by which the hydrogen donating water molecules move across the first and second shells, thereby causing this fluctuationopen